Observational constraints on the linear fluctuation growth rate 
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Several experiments in the near future will test dark energy through its effects on the linear 
growth of matter perturbations. It is therefore important to find simple and at the same time 
general parametrizations of the linear growth rate. We show that a simple fitting formula that 
generalizes previous expressions reproduces the growth function in models that allow for a growth 
faster than standard, as for instance in scalar-tensor models. We use data from galaxy and Lyman-a 
power spectra to constrain the linear growth rate. We find 7 = 0.61q 3 for the growth rate index 
and 77 = O.OIq j for the additional growth parameter we introduce. 
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INTRODUCTION 



After several years from the first works [l|, the evidence of dark energy (DE) still rests primarily on background 
quantities like the luminosity distance and the angular diameter distance 0, [j, [B| . Only recently the cross-correlation 
of the integrated Sachs- Wolfe effect with the large scale structure yielded an independent proof of the existence of 
£V. dark energy that rely on the linear growth of the gravitational potential [fj]. In the future, galaxy and Lyman-a 
power spectra at high redshift and weak lensing surveys from ground and from space will offer the opportunity to test 
competing dark energy models to a very high precision using a mix of background, linear and non-linear indicators 

In order to use the growth rate of linear fluctuations as a test of DE, it is however necessary to formulate simple 
expressions that embody a large class of models, so as to provide observers with a practical tool to analyse the data. 
This procedure has proved most convenient with the equation of state of DE, whose parametrizations afford the 
community to quickly compare the different experiments and to optimize the design of further surveys (see eg [§]). 

Since many years it has been known that a good approximation to the growth rate of the linear matter density 
contrast 6 in standard models of gravity and dark matter is given by the simple expression [TcL fll] | 
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= j^^L, (i) 



where the prime stands for derivative with respect to log a and where 7 ~ 0.55. In [12| it was found that this growth rate 
also accounts for coupled dark energy at small z and small coupling (3 (see definition below) with 7 ~ 0.56(1 — 2.55/3 2 ). 
The formula works quite well also in the case in which DE is described by a mildly varying EOS w{z) when generalized 
with [13 7 « 0.55 + 0.05(1 + w(z = l)). 

However, this formula has an obvious drawback: assuming a constant or weakly varying 7 > 0, it implies s < 1 
5_j 1 at all epochs since in standard cosmology the matter density parameter is always < Sl m < 1. More exactly, if 7 
is a constant, then s is always either smaller or larger than unity. Therefore, this fit is unable to test for deviations 
from the standard paradigm of fluctuation growth, which assumes s = 1 for the matter-dominated universe at high 
redshifts. This assumption is of course dangerous in view of the multitude of dark energy or modified gravity models 
currently being studied and should not be taken for granted without serious scrutiny. 

In this paper we show first that the s > 1 behavior takes place in one of the simplest class of modified gravity 
theories, namely scalar-tensor models or their Einstein frame counterparts, and then we introduce a simple new fitting 
formula that generalizes eq. ([I]) . We note that a similar faster-than-standard growth has been found also in TeVeS 
modified gravity models [LJ]. Finally, we compare our fitting formula with the (scanty) data available at the present. 



II. THE GROWTH FUNCTION IN SCALAR-TENSOR THEORIES 



In this section we show that generically s > 1 at early times and s < 1 at late times in a simple class of modified 
gravity theories. 

Let us consider a generic scalar-tensor model in Einstein frame. It is well-known that in the Einstein frame, the 
frame in which the gravitational equations are in the Einsteinian form while matter is not conserved, the scalar degree 
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of freedom of a scalar-tensor theory acts as a new force on matter, as expressed by the new conservation equations 
(here we use subscripts c for CDM, b for baryons, <jf> for the scalar field, and 7 for radiation ) [TR HiL Il7j 

T "Mc) = -Cc(tf)T (c) *. M , (2) 

T u;n(b) = - C b{<t>)T(b)<l>;», (3) 

T u-M4>) = l°bWb) + C C {4>)T {C) ]^ , (4) 
T M , > = 0, (5) 

where for generality we assumed that the coupling functions depend on the species (while radiation remains uncoupled 
because it is conformally invariant). In a FRW metric with scale factor a these equations become (we assume flat 
space throughout) 



cp + ZH^ + Vj, = k{(3cPc + PbPb) , 

p c + 3Hp c = -K/3 c p c <j) , 

p b + 3Hp b = -Kf3 b pb4>, 

p'-y + 4Hp~ ( = , 

3H 2 = K 2 (p b + p c + P4,), (6) 

where k 2 — 8nG, (3 C — C c /n, f3 b — C b /n, H = a/a (note that we use a coupling (3 which is \/2/3 the (3 used in ref. 
[Trip. Then we have immediately 

Pc.b = p{0)c, b a~ 3 exp I - J (3 c , b {4>)d<j) \ . (7) 

To simplify the analysis and to satisfy local grav ity c onstraints we put from now on fib — and /3 C — (3 — const. 

As it has been shown in refs. [T3|, for (3 < y/3/2 the standard matter era that precedes the final acceleration is 
replaced in this coupled model by an epoch in which the energy density fractions fl m ,il,^ of matter and field are 
constant and equal to 

^0 = \P 2 , (8) 

and fl m = l — fi^. During this epoch one has 4>' — 2/3 (the prime stands for d/dloga) and the scale factor grows 

as a ~ 1 3 (!+™=) with w e — 2/3 2 /3 (these values are approximated since are obtained neglecting both baryons and 
radiation). This new matter era has been denoted as 0MDE in previous works. This occurs when the potential is 
negligible with respect to the field kinetic energy. Since the potential is dominating the final accelerated epoch, it is 
clear that the (/)MDE generically will take place before acceleration and, of course, after the radiation era. This in fact 
is what has been observed in several numerical and analytical investigations, for instance in the case of exponentials 
and inverse power-law potentials V{4>) — A$~ n [li§]. In the Jordan frame, where matter is conserved, the </>MDE 
corresponds to the standard solution of the Brans-Dicke original theory (which is derived in absence of a potential) 



upon the substitution 



aj ~ , (9) 



Therefore, the <£MDE is quite a generic feature of scalar-tensor models and it also shows up in some f(R) models 
We now show that during </>MDE the growth of fluctuations is faster that in a standard matter era. In ref. [131 
it has been shown that the perturbation equations in the sub-horizon regime is 

1 + ji - /^') 6 'c - \{icc^ c + i bc s b n b ) = , (11) 

where again the prime stands for derivation with respect to log a, and where jij = 1 + 2(3i[3j and H is the conformal 
Hubble function TL = aH . We assume the baryon component to be negligible. Then the fluctuation equation can be 
solved analytically during the </>MDE: 

S ~ a 1+2 P 2 , (12) 
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from which it appears that s = 1 + 2/3 . In the Jordan frame, the growth rate is instead 

1 + 2/3 2 

8 j ~ a"^ ~ a 1 ^ , (13) 

which also gives a rate larger than unity. When the </>MDE ends and acceleration takes over, the rate s declines 
steadily to zero as in standard cases. Therefore, as we anticipated, s goes from a value larger than unity to a value 
smaller than unity. 

III. A GENERALIZED FIT 

We now proceed to find a convenient fit to the full evolution of 6(a) for the coupled models introduced in the 
previous section. In the standard scenario 6 obeys the equation 

6"(a) + (1 + j{)6'(a) - ~O ro £(a) = 0, (14) 



where 



S = 4(1 + 3ttf,(a)fi,(a)) = - J(l + W 11 ) (15) 



The solution can then be approximated as 



In our modified gravity eq. lfT4|) becomes 



£( a ) =e /o°&«'nm{«T - (16) 



5 " + ( 1 + % ~ PA 5 ' ~ i"m(l + 2/? 2 )<J = ■ (17) 



n ) 2 

One simple possibility would be to generalize lfl6|) as 

S( i a)=eSo d «' r ^W , ( i +^\ (18) 

with c a parameter to be determined by a least square fit. The choice of a (3 2 behavior is suggested by the fact that the 
0MDE depends only on (3 2 . However, this new fit is not very practical because it contains the function fi m (a) that 
should be obtained by numerically integrating the background equations and therefore depends on the field potential. 
To overcome this difficulty, we propose to use instead the standard expression for fi m : 

= o z <?" o ^ -*» > ( 1Q ) 

" m ,o + (1 - " m ,o)a J! " 

where w = (loga) -1 J w(a)da/a and the subscript denotes the present time. For the coupled dark energy model we 
are considering here, we approximate w(z) ss w^(z = 0); although one could easily expand w(z) to higher orders, our 
approximation is sufficient to show that our generalized fit works well. Therefore we define the rate 

S ^n£(ay(l + c(3 2 ). (20) 

In this way, the growth rate can be parametrized by fi mj o,7, and the combination r\ = c/3, plus the parameters that 
enter w(z). Now, even in the limit fi m — > 1 one has s / 1. In the next section we show that this generalized fit is 
indeed a good approximation. Since we know that during the </>MDE (i.e. at high z, for which 0,$ w 1) one has 
s = 1 + 2/3 2 we can anticipate that the result will be close to c « 2. 



Concluding this section we note that eq. I|20|) should be seen for what it is, i.e. a phenomenological fit. The 
relation of Ct m ,o, w(z), 7 and 77 to the underlying theory will of course depend on the theory itself. For instance, 
the identification of £l m ,o with the presently clustered mass in galaxies and clusters of galaxies is actually a model- 
dependent assumption; if gravity is not standard this assumption is likely to be incorrect. All we are assuming here is 
that the Friedmann equation can be written as the sum of two components, one that dilutes as £l m ,oa~ 3 and the other 
as (1 — £! m ,o)a~ 3 ' 1+tu ' ) ; if the gravitational equations are not standard, one has to define £l m (a) such that the above 
parametrization is still valid. The advantage of using (|19H20p is that both background and linear growth are fitted 
by the same expression for w(a); that is, once one adopts a prescription for w(a) one can fit all the data by simply 
adding the two parameters 7 and r\ (plus possibly further parameters to account for the anisotropic stress, see eg. 
[2l|). Of course, in principle one could proceed in different ways: for instance, one could parametrize fi m (a) so that 
values larger than unity in the past were allowed so as to force s > 1. Trivially, in fact, our parametrization above 
could be written equivalently defining a new density TO = ttm\l + ry) 1 ^ 7 ; this density parameter would however not 
be the same that appears in the Friedmann equation. 
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Figure 1: We compare the functions gfu = Sfu/a (red solid curves), given by the fitting formula (18)1 with the best fit 
(parameters 7=0.56, c = 2.1) for two different values for /3, to the exact solutions gexact = S exa ct/a (black dashed curves) of the 
differential equation lfl4|| for the growth rate. The curve for (3 — also gives the standard best fit (i.e. for c = 0). All curves 
are normalized at unity at z = 20. The cases f3 = 0, 0.2 are for an exponential potential, the case /3 = 0.1 for V ~ l/4>- 



IV. COMPARING THE FIT TO THE NUMERICAL RESULTS. 

We solved numerically the background equations of the system |(6|) choosing an exponential form for the potential 
V{4>) and neglecting the fraction of baryons and radiation. Then we solved numerically the differential equation (fTTl) . 
thus obtaining a solution (that we denote 5 exac t) which depends on the value of the coupling constant (3, For (3 ranging 
between and 1/2, we found that the values of the parameters 7, c appearing in (|20|) . which produce the best fit to 
Sexact are 7 = 0.56 and c = 2.1. Our best fit is therefore 

^(a,/3)^e/o Q ^'^ , («)°- 56 (i+2-i/5 2 ) j (21) 

where we remark again that we use the standard expression for fl m (a). 

This new function is really a good approximation to the exact solution S exact , for different values of the coupling 
constant (3, as one can see in FigQ] where the curves of the growth factor g = S/a for two different (3 and the 
corresponding gju are plotted. In Fig. 2 we present the level of accuracy of the fitting formula. As it can be seen we 
find fits to better than s=sl% for different values of f3. Moreover, we find that the best fit values of the parameters do 
not depend on the actual value of the present matter density O m ,o- We experimented also with an inverse power-law 
potential and found that also in this case eq. (|2T]l is a good fit (see curve for (3 = 0.1 in Fig. Q]). Without the 
^-correction the relative error (Sfu — & exact) I & exact becomes larger than 15% already for (3 — 0.2. 

V. COMPARING THE FIT TO THE OBSERVATIONS 

In the previous section we have seen that the expression 

s = n$(ay(l + r,), (22) 

where £1$ (a) is given by eq. I|19|) gives a good fit to the evolution of S during both the decelerated and accelerated 
regimes for coupled dark energy models with constant (3 < 1/2 if 7 w 0.56 and 77 = 2.1/3 2 . Here we take some 
preliminary steps towards comparing the fit f22|) to the observations. An indication for a positive rj could signal an 
attractive force additional to standard gravity as in a scalar-tensor model; on the other hand one can speculate that 
a negative 77 could be related to the slowed growth induced by a hot matter component. 

We consider the following data: a) Lyman-a power spectra at an average redshift z = 2.125, z = 2.72 [2a ] . z = 3 
[2~3 |: b) the normalization as inferred from Lyman-a at z ranging between 2 and 3.8 [25]; c) galaxy power spectra at 
low z from SDSS [2§| and 2dF (2?J • From the three Lyman-a and the SDSS spectra we estimate the ratios 
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Figure 2: Level of accuracy of the best fit to the exact solution for the growth rate. For /3 ranging between and 0.3 the fits 
are better than 1.2%. Without the correction the relative errors would be larger than 15% already for /3 = 0.2. 



for the values of ki for which there are tabulated value of the spectra (or for interpolated values and errors when the 
tabulated wavenumbers differ). For the as data we estimate the ratios between successive values of z, 

r(zi ' z2) = 5li (24) 

<7 8\ z 2) 

(note that ref. [25j reports the values of as extrapolated at the present epoch ). 

For the Lyman spectrum at z — 3 and for 2dF (z = 0.15), the authors of [13, H3] give directly their estimation of the 
growth rate, s b s =0.49±0.10 for 2dF and s b s =1.46±0.29 for the Lyman-cv data. Then we compare the observations 
to our fit by using the likelihood function 



JVexp) f--^ TTjT ex P>. -— KIT . (25) 




where the errors a% are obtained from the quoted errors on P(k) and as by standard error propagation. 

As we will see the data available at the present are not sufficient to set stringent limits to the growth function. 
Moreover, there are several sources of possible systematic effects that we cannot account for. For instance, the matter 
spectra derived from Lyman-cv clouds are obtained through calibration (ie. bias correction) with iV-body simulations; 
these simulations have been generated only for a limited set of cosmological models. It is difficult to quantify the 
impact of this limitation upon our results; the fact that we consider ratios of spectra from similar sources (eg Lyman-a: 
clouds) might however alleviate the problem since one can expect that the calibration errors are only weakly dependent 
on redshift. For this reason we consider separately the ratios of the high-z Lyman-a spectra to the low— z SDSS galaxy 
spectra; our final results do not take these into account. 

The current observational situation is summarized in Fig. J3]) (and the associated Table I), in which we plot the 
data we used in this work, along with the ACDM growth rate and with our best fit (see below). This figure gives a 
clear idea of the potential for improvement in the observational estimation of the growth rate. 

We assume that the function s depends on four parameters, (f2 m ,o, wo, 7, 77). We assume also a flat prior O m ,o 6 
(0.05,0.4) and wq e ( — 1,-0.6) which generously accounts for the supernovae constraints (neglecting the phantom 
region). Our main result is contained in Fig. |(5]), which displays the likelihood contour plots at 68%, 95 and 99.7% 
in the plane (7,77), marginalizing over O, m ,o,wo. Remarkably, the best fit values practically coincide with the ACDM 
prediction, (77,7) = (0,0.6). However the likelihood extends considerably on both negative and positive 77 and even 
negative values of 7 are not excluded beyond 3a. In Figs. ((MZj) we plot the marginalized ID likelihoods for 7 and rj. 
The results are tabulated in Table II. The best fit values and la errors are 

7=0.6018$, 77 = 0.00^. (26) 

As we anticipated, the current data impose only very weak constraints on 7,77 . For completeness, we also quote in 
Table II the best fit and errors on ^standard, i.e. assuming a standard model in which 77 = 0. Even in this case the 
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z 


s 


ref. [23] 


2.125; 2.72 


0.50; 0.98 


ref. [25] 


2.2; 3 


-1.147; 1.175 


2.4; 3.2 


-0.94; 1.198 


2.6; 3.4 


-0.686; 2.010 


2.8; 3.6 


-0.908; 1.778 


3; 3.8 


-1.207; 1.799 


ref. [24] 


3 


1.46±0.49 


ref. [27] 


0.15 


0.49±0.10 



Table I: Summary of observational data. We report in the z and s columns either the corresponding ranges or the central value 
and errors. For the as data or ref. f 273) we chose to report the errorboxes on s obtained using the ratios at the given redshifts. 
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2(7 
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V 


o.ootg:3 


+0.58 
-0.38 


+1.1 

-0.58 


7 


o.6o±r 30 


+0.97 
-0.49 


+ 1.6 
-0.74 


^/standard 


0.60t°;« 


+0.77 
-0.40 


+ 1.4 
-0.50 



Table II: Best fit and errors (marginalized over all other parameters). 

likelihood distribution for 7 remains very broad, although now negative values are rejected at more than 3a. Including 
the ratio of Lyman-a to SDSS power spectra has a minor effect on 7 and moves the best fit of 77 to —0.2. 
Assuming 77 < 0.58 at 2a we can derive an upper limit to the coupling (3 introduced in Sect. 2, 

[3 < 0.52 (27) 

(at 95% c.L). This limit is very weak when compared to the CMB limits [l8| but it is nevertheless interesting since 
it is independent and derived uniquely from the growth rate at small redshifts. 

VI. CONCLUSIONS 

The search for useful parametrizations of the dark dynamics is important since as it has been shown several times 
every parametrization introduces some arbitrariness in the way data are analysed |29j|. In particular, with the advent 
of models of dark energy based on modification of Einstein's gravity, we have become aware of many possible trends, 
both at the background and at the perturbation level, that are not easily accounted for with earlier parametrizations. 
In this paper we introduced a generalized form of parametrization of the growth rate that allows for a rate s^l, i.e. 
faster or slower than the standard matter-dominated growth. We show that this parametrization is suitable to model 
the fluctuation growth in coupled dark energy models and in scalar-tensor models. 

We have analysed the current data in search of observational constraints on 7, 77. Considering data from Lyman- 
a and galaxy power spectra at various redshifts we have obtained (rather weak) constraints on both parameters. 
The best fit turns out to be very close to the ACDM predictions. Many future experiments based on weak lensing 
or baryon oscillations will be able to estimate the growth rate and other fluctuation parameters with much higher 
precision [2]], [H, [28| . We expect therefore that the constraints derived in this paper will soon be superseded by much 
more precise ones and that new estimates of the growth factor will help in clarifying the nature of dark energy. 
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Figure 3: Summary of experimental data for the growth rate s, as detailed in Table I. The big coloured errorboxes represents 
the ratios a(zi)/cr(z2) for various z intervals of ref. [25j] (three additional very large errorboxes have been excluded from the 
plot but not from the analysis); the smaller black box represents the average spectral ratio for the Lyman-a data of ref. [231 ]. 
The two points with errorbars are from ref. [27| and ref. [24j]. The black solid line is the ACDM model, the red curve is 
the coupled dark energy model with = 0.2 and /3 = 0.4 (i.e. r\ = 0.34) and the dashed blue curve is the overall best fit 
(fim.o, wo, 7, 17) = (0.05, -0.6, 0.4, 0.45). 




Figure 4: Tridimensional likelihood function marginalized on fl m ,o and wo ■ The peak corresponds to (7, 7?) = (0.6,0). 
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